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We present a discussion of the effects induced by the bulk viscosity on the very early Universe 
stability. The matter filling the cosmological (isotropic and homogeneous) background is described 
by a viscous fluid having an ultrarelativistic equation of state and whose viscosity coefficient is 
related to the energy density via a power-law of the form £ = Cop" ■ The analytic expression of the 
density contrast (obtained for v — 1/2) shows that, for small values of the constant (o, its behavior 
is not significantly different from the non-viscous one derived by E.M. Lifshitz. But as soon as 
£o overcomes a critical value, the growth of the density contrast is suppressed forward in time by 
the viscosity and the stability of the Universe is favored in the expanding picture. On the other 
hand, in such a regime, the asymptotic approach to the initial singularity (taken at t = 0) is deeply 
modified by the apparency of significant viscosity in the primordial thermal bath i.e. the isotropic 
and homogeneous Universe admits an unstable collapsing picture. In our model this feature regards 
also scalar perturbations while in the non-viscous case it appears only for tensor modes. 

PACS numbers: 98.80.-k, 51.20.+d 

I. INTRODUCTION 

The very early stages of the (homogeneous and isotropic) Universe evolution, say at tem- 
peratures greater than O(10 16 GeV), are characterized by a thermal history which can not 
be regarded as settled down into the equilibrium. Indeed at sufficiently high temperatures 
the cross sections of the microphysical processes, responsible for the thermal equilibrium, 
decay like 0(1/T 2 ) and they are no longer able to restore the equilibrium while the Uni- 
verse expands. Thus, going backward, we meet stages where the expansion has an increasing 
rate and induces non-equilibrium phenomena in the primordial bath. The average effect of 
having a microphysics unable to follow the Universe expansion by equilibrium stages, re- 
sults into dissipative processes appropriately described by the presence of bulk (or second) 
viscosity. 

As shown in [3, [3, Q the thermodynamical consequences which bulk viscosity induces on 
the very early Universe history, can profoundly modify its dynamics and configurations in 
which matter is created near the singularity become possible. Such scenarios arise because 
the bulk viscosity coefficient £ is a function of the Universe energy density p and can be 
modeled in terms of a power-law as Q = Ow 1 ' where Co, v = const (for a discussion of this 
constitution equation see [1] and see also jf| 

Aim of the present work is to investigate the effects that bulk viscosity has on the stability 
of the isotropic Universe, i.e. the dynamics of cosmological perturbations is analyzed when 
viscous phenomena affect the zero and first order evolution of the system. We consider a 
background corresponding to a Robertson- Walker model filled with ultrarelativistic viscous 
matter, whose coefficient £ corresponds to the choice v = 1/2 and then we develop a 
perturbation theory which generalizes the Lifshitz works 0, Q to the presence of bulk 
viscosity. Though the analysis is performed for the case of a flat model, nevertheless it 
holds in general as soon as the perturbations scale remains much smaller than the Universe 
radius of curvature. In this respect we deal with perturbations such that rjq <C 1, ^n/q being 
the size of the coordinate scale and r\ the conformal time variable. Since the dynamics we 
consider holds near the singularity for r\ <C 1, then we make allowance for arbitrarily large 

1 /2 

values of q and therefore the condition for the general validity q ^> 2tt\1C\ (/C being the 
Robertson- Walker curvature parameter) can be always fulfilled. 

As issue of our analysis we find that two different dynamical regimes appear when viscosity 



2 



is taken into account and the transition from one regime to the other one takes place when 
the parameter £ overcomes a given threshold value. However in both these stages of 
evolution the Universe results to be stable as it expands; the effect of increasing viscosity is 
that the density contrast begins to decrease with increasing n when fo is over the threshold. 
It follows that a real new feature arises with respect to the Lifshitz analysis when the 
collapsing point of view is addressed. In fact, as far as £0 remains below the threshold value, 
the isotropic Universe approaches the initial Big-Bang with vanishing density contrast and 
its stability is preserved in close analogy to the non- viscous behavior. But if £0 overcomes its 
critical value, then the density contrast explodes asymptotically (77 — > 0) and the isotropic 
Universe results unstable approaching the initial singularity. In the Lifshitz non-viscous 
analysis this same backward in time instability takes place only when tensor perturbations 
(gravitational waves) are taken into account, since their amplitude increases backward as 
the inverse of the cosmic scale factor. 

The new feature induced by the bulk viscosity consists of having instability simply in 
correspondence to scalar perturbations induced by fluctuations in the matter filling the 
Universe. The cosmological interest in such instability of the primordial Universe (toward 
scalar perturbations) comes out reversing the picture from collapse to expansion and taking 
into account the time reversibility of the Einstein equations. In fact if the early Universe 
does not emerge from the Planck era peaked around the Robertson- Walker geometry (indeed 
a good degree of generality in its structure is predicted either by classical and quantum 
argumentation 0]) then it can not reach (according to our analysis) an homogeneous and 
isotropic stage of evolution before the viscous effect become sufficiently small. 

Since a reliable estimation [lOj fixes the appearance of thermal bath into the equilibrium 
below temperatures O(10 16 GeU) and this limit corresponds to the pre-inflationary age, our 
result supports the idea that an isotropic universe outcomes only after a vacuum phase 
transition settled down. In Sec. II we present the paradigm underlying perturbations 
theory to the Robertson- Walker Universe and provide the key expressions to determine 
the inhomogeneous corrections corresponding to different states of the matter filling the 
cosmological space. In Sec. Ill a brief review of the Lifshitz analysis is presented in order 
to compare our subsequent results with those ones outcoming from a non-viscous model. 
In Sec. IV we provide the basic equations governing the perturbations dynamics in a 
viscous isotropic Universe. In particular we generalize the energy-momentum tensor of the 
cosmological fluid in view of including the viscous effects and then derive the new zero and 
first order dynamics; the Friedmann equations corresponding to the viscous case are solved 
for the choice v = 1/2. In Sec. V we solve the perturbations equations for the same value 
of the parameter v\ the behavior of the 3-metric perturbations and of the density contrast 
is provided in the asymptotic limit to the Big-Bang and the outcoming issues are discussed. 
In Sec. VI brief concluding remarks follow and possible upgradings for the description of 
the bulk viscosity cosmological effects are discussed. 



II. PERTURBATIONS THEORY TO THE EINSTEIN EQUATIONS 

In order to describe the temporal evolution of the energy density small fluctuation, we 
develop a perturbations theory on the Einstein equations. We limit our work to the study 
of space regions having small dimensions compared with the scale factor of the Universe 
a [111 ]. According to this approximation, we can consider a 3-dimensional Euclidean (time 
dependent) metric as spatial component of the background line element 

ds 2 = dt 2 - a 2 (dx 2 + dy 2 + dz 2 ) . (1) 

In the linear approximation, perturbed Einstein equations are 

5R» - l -5^5R = ZwGSTZ , (2) 

where the term 6T£ represents the perturbation of the energy-momentum tensor which de- 
scribes the properties of the matter involved in the cosmological collapse. The perturbations 
of the Ricci tensor SR^ can be derived from the metric perturbations = —Sg 1 ^, since the 
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general expression for the perturbed curvature tensor is 

dRfwp = 2 (^V; p; v hp- p; y — ^pp; 17 — ^p; f; p — ^"j p; p fytw; p) ' (^) 

For convenience let us introduce a new temporal variable ry, given by the relation dt = a dr), 
and use the symbol (') for its derivatives; we moreover impose, without loss of generality, 
that the synchronous reference system is still preserved under perturbations, so that 

/ioo = h 0a = . (4) 

If we consider the background metric |T]) the perturbations of the mist components of the 
Ricci tensor and of the curvature scalar read: 

= + K\l ~ m ~ hi) - ^ hi" -^hZ-± h'SP a , (5b) 

5R =-l(hl^-h^)-^h"-^h' , (5c) 

By using these expressions we are able to rewrite the left-hand side of Einstein equations 
through the metric perturbations hp* in order to develop the perturbations theory after 
describing the matter properties via an appropriate energy-momentum tensor. 



III. PRIMORDIAL UNIVERSE IN ABSENCE OF VISCOSITY 



In this Section we assume that the Universe, in its primordial expansion, behaves like a 
perfect fluid. This hypothesis was made by Lifshitz in his works @, H[ and can be expressed 
writing the energy-momentum tensor in the form: 

= {p + p)u tl u u - pg/ , (6) 

where p, p are respectively the pressure a the energy density of the fluid, and is its 
4- velocity expressed in the comoving system we consider, i.e. 

u° = l/a u a = 0. (7) 

Using the synchronous character of the perturbed metric, we are now able to write the 
perturbations of the above energy-momentum tensor: 

6T° = Sp, ST a = a(p + p)Su a , ST* = v 2 s Sp , (8) 



v s = \J Sp/Sp being the sound speed of the fluid. 

Since we use an Euclidean background metric, we can expand the perturbations in plane 
waves of the form e tq r , where q is the adimensional comoving wave vector being the physical 
one k = q/a. Here we investigate the gravitational stability which is described by the 
behavior of the energy density perturbation expressible only by a scalar function; thus we 
have to choose the scalar representation of the metric perturbations in order to involve a 
change, not only in the gravitational field, but also in the velocity and in the energy density 
[1, Oj]]. Such a development is made by the scalar harmonics Q — e lq r , from which the 
following tensor 



Qi - \€Q , H 



Q ■ (9) 



can be constructed. We can now express the time dependence of the gravitational pertur- 
bations through two functions A(?y), p(r)) and write the tensor h@ in the form 

hi = Hv)P£ + ^l)Qi , h = p( V )Q . (10) 
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Let us now consider the primordial stages of the Universe expansion, i.e. 77 -C 1, when 
the radiation-like density dominates the matter one. The equation of state is p = p/3, from 
which the relations (for a flat Universe JC = 0) arise 



p = Ca 4 , a = a\r\ , v 2 s — 1/3 , (11) 



where C is an integration constant and a\ = W8irGC/3. In this approximation we can 
obtain the basic equations which describe the temporal evolution of the perturbations. 
Expressing equations {5| through the representation (jTDJ) and using expressions ([5]) in the 
form 8T@ = — 5@v 2 STq, the perturbed Einstein equations give, respectively for a ^ (3 and 
for contraction over these indexes, two equations for the metric perturbations 

A" + -V-^(A + /i) = ! M " + ^ M ' + ^!(A + M ) = 0. (12) 

r\ 3 r\ 3 

Furthermore, taking the 0-0 components of ([2|), we can express the energy density directly 
from the adopted functions A and p in the form 



r Q 
dp 



24nGa 2 



q 2 (\ + p) + —p' 
a 



(13) 



Among the solutions there are some which can be removed by a simple transformation 
of the reference system (compatible with its synchronous character), and therefore they do 
not represent any real physical change in the metric. The corresponding expression for the 
metric perturbations can be established, a priori, through a coordinates transformation [ll[ 
taking into account the constraint (@|: 

K = fh%j d ^ + ^h€ + {fk p + f p a ) , (i4) 

where /o, f a are arbitrary (small) functions of the coordinates. 

In the assumption r\q <C 1 the equations (|12|) admit, in the leading order, the solutions 

A=— + C 2 , t i = -^-C 1 r 1 +C 2 , (15) 
77 3 

where the fictitious solutions fj 14[) . which in our ultrarelativistic approach assume the form 
A — /j = const (/o = 0, f a = P a ) and A + p ~ 1/rj 2 (/o = Q, f a = 0), are excluded. The 
final expressions for the gravitational perturbations and for the density contrast Sp/ p can 
be obtained substituting this solutions in (p~0|) and (fT3|) 



- = ^(Civ + C 2V 2 )Q. (17) 
p 9 

Here the constants C\ , C2 must satisfy the conditions expressing the smallness of the per- 
turbations at the moment 770 when they arise; assuming that harmonics Q are of the unity 
order magnitude, the inequalities A <§; 1, fj, <C 1 give the constraints C\ <§; t\q <C 1 and 
C 2 <1. 

The expression of the cosmological perturbation (fT?} contains terms which increase, in 
an expanding Universe, proportionally to positive powers of the scale-factor a = air]. This 
expansion can't, nevertheless, imply the gravitational instability: if we consider the magni- 
tude order -q ~ the conditions satisfied by the constants Ci, C2 imply that the density 
perturbation remains small even in the higher order of approximation. This behavior of 
the cosmological fluctuation yields the gravitational stability of the primordial Universe; 
the only stability we can found in a non-viscous Universe [8[ is provided by the tensor 
perturbations /i^and takes place approaching backward the Big- Bang. 
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IV. DYNAMICAL DESCRIPTION OF THE VISCOUS UNIVERSE 



In the last Section we assumed that the Universe could be described by the energy- 
momentum tensor of an ultrarelativistic perfect fluid. The immediate generalization is 
to consider the presence of dissipative processes within the fluid dynamics as expected at 
temperatures above O(10 16 GeU); this correction is represented by an additional term in the 
expression of the energy-momentum tensor ^ and it can be derived from thermodynamical 
properties of the fluid (in particular from the law of increasing entropy (T3, Il3l|). Using the 
conservation law T" v = we arrive [3, HH to express the energy-momentum tensor of an 
ultrarelativistic viscous fluid in the form 



T, 



(p + p)Uf,u u ~V9p,v 



P=P- 



(18) 



where p denotes the usual thermostatic pressure and £ is the bulk viscosity coefficient. In 
this work we neglect the so called shear viscosity (first viscosity) since in the case of isotropic 
cosmological evolution there is no displacement of the matter layers with respect to each 
other and this kind of viscosity represents the energy dissipation due to this effect. 

The coefficient £ is not constant and we have to express its dependence on the state 
parameters of the fluid. In the homogeneous models this quantity depends only on time, 
and therefore we may consider it as a function of the Universe energy density p. According 
to literature developments [l|, 0, H, H[ we assume that £ depends on p via a power-law of 
the form 



C = Co/ 



(19) 



where Co is a constant and v is an adimensional parameter. The behavior of this parameter 
is derived by V.A. Belinskii et al. [lj for asymptotic values of the density energy yielding 
the constraint 0^^^1/2in the region of large p. 

Let us now perturb the viscous energy-momentum tensor in our synchronous reference 
system obtaining the expressions 



5T°=5p, 5T«=a{p + p)5u a , ST? = 5* 



-C 2 s Sp + C[Su\ + —r 



ft' 

2a 2 



and the relation between the mist components of the tensor 

b! 



ST? = 



2a a(ft + p) 



(20) 



(21) 



using the background metric |T]) and the expressions ([7|) of the 4- velocity u M in a comoving 
system, the viscous pressure p becomes 



P = P- 



3C^ 



(22) 



and the quantity c s is given by 



3C0 -n VP 



(23) 



Of course the presence of viscosity does not influence the expression of the Ricci tensor 
and its perturbations, thus we can still keep the expressions (0 and use the formula (f2Tj) to 
build up the equations which describe the dynamics of h^ a and Sp. It is convenient to choose, 
as final equations, the one obtained from the Einstein ones for a ^ f3 and for contraction 
over a and /3, which read respectively 



2a' 



(hl? J + h%Z-h&-h£)+hZ" + — ^'=0, a^p, 



(24) 



-(h%°-h%)(l + 3<*) + h"+ 



-(2 

a v 



3c 2 - 12ttG- C) ti 
a' ' 

3C , L . a < 



2a(p + p) 



(25) 



hi, a 
"0,7 



G 



Furthermore the fictitious solutions (| 14[) stand also in presence of dissipative processes 
because they are founded by a transformation of synchronous reference system. 

Since we want to describe the gravitational instability in presence of viscosity, we consider 
the scalar representation (flT)]) which, once substituted in the last expressions, yields the 
equations describing the perturbations temporal evolution. Let us now express the time 
dependence of the model variables; as in the last Section, we consider the earlier stages of a 
fiat Universe corresponding to 77 <C 1 and with the equation of state p = p/3. The Universe 
zero-order dynamics is described by the equation of energy conservation and the Friedmann 
one, which are respectively 



p' + 3-(p + p)=0, - = J-nGp. (26) 
a a z V o 

In order to integrate these equations we assume, according to the large energy density of 
the primordial expansion, v = 1/2 [j| Substituting (f22|) into the above equations we 
obtain, for v = 1/2: 



p = Ca- {2+2u \ a^am 1 ^, w = 1 - V54irG Co , (27) 

being C an integration constant and a\ = (8w 2 7rCG/3) 1 / 2 ". Since we consider an expanding 
Universe, the factor a must increase with positive power of the temporal variable thus we 
obtain the constraint < u> ^ 1. 

Using these explicit dependences we get two equations for the A, p time functions 

A" + — A'-£ (A + M )=0, (28) 

LOT] 3 



q 2 Co 



WC/3af - 3Co/w 



(29) 



whose solutions describe the evolution of the gravitational perturbations and the energy 
density fluctuation. 



V. THE BEHAVIOR OF DENSITY PERTURBATION 

In the model developed in this paper, we study the gravitational collapse dynamics of 
the primordial Universe near the initial Big-Bang for 77 <C 1. As in Lifshitz's works, we now 
analyze the case of perturbations scale sufficiently large to use the approximation rjq <c 1. 
In the non-viscous model the cosmological stability of the isotropic Universe is guarantied 
by the positive power-law exponents of the density contrast evolution and by the constraint 
for the constants C\, C^. In our scheme equations (|2"5]) and (|29[) admit asymptotic analytic 
solutions for the functions A and p; in the leading order A takes the form 

X = ^T + C ^ (3°) 

where Cj., Ci are two integration constants. Substituting this expression in (|29[) we get, in 
the same order of approximation, the behavior of the function p as 

^=^3+^: (31) 

where we have excluded the non-physical solutions (|14p as written in the form \ — p = const. 
The constant C\ is given by the expression C\ = A/B(3 — l/u>), A and B being constants 
involved in the equation for p having the form 

A + ^-W^ , B = 12 -f 5" . (32) 

*1 



4VC/3< - 3Co/^ a 
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Let us now write the final form of the perturbations pointing out their temporal depen- 
dence in the viscous Universe. The gravitational perturbations (jTUJ) become 

K = ^tf + ^Ql + C2(QZ + rS) ■ (33) 



and the density contrast reads 

'C 1V 3 - 2 ^ + C 2 v 2 + C 3V 3 - 1/u + , (34) 



5 -P-F 

— 1 UJ 



where C3 = 3A/q 2 coB, and F u = uj 2 Qq 2 /9. As in the non-viscous case, we now impose the 
conditions expressing the smallness of perturbations at the initial time 770 ■ The inequalities 
ftfcl and Sp/p <gC 1 yield only two fundamental constraints for the integration constants: 
Ci < tfj^ 1 and C 2 < 1 for any lo- value within the interval (0, 1]. Furthermore we find an 
additional condition which involves the wave number q and the integration constant C; in 
particular a rough estimate for lo < 1/3 of the inequalities C\ <C Vo^~ 3 an d C3 <C r^/ u 3 
yields the condition q <C (GCr/o) -1 ' 2 " which ensures the smallness of the cosmological 
perturbations. 

Using the hypothesis r\ <C 1 we can get the asymptotic form the corrections to the 
cosmological background. The exponents of the variable 77 can be positive or negative 
according to the value of the viscous parameter lo, which is always positive but less than 
unity. This behavior produces two different regimes of the density contrast evolution: in 
the first case 2/3 < lo ^ 1, the perturbation increases forward in time which corresponds 
qualitatively to the same picture of the non- viscous Universe where the density fluctuation 
increases but remains small; while in the other case < lo < 2/3 the density contrast is 
suppressed behaving like a negative power of 77. When the density contrast results to be 
increasing, the presence of viscosity induces a damping of the perturbation evolution in 
the direction of the expanding Universe, so the cosmological stability is fortified since the 
leading 77 powers are smaller than the non- viscous ones (|17[) . 

For < lo < 2/3, i.e. for a large coefficient £o; the density contrast evolution is deeply 
modify by the presence of viscosity. We get the leading order expression 

^ = , k = 2/u> - 3 > . (35) 

This behavior is very different from the non- viscous one since it yields a strong damping of 
the density contrast. In this regime the density fluctuation decreases forward in time but 
the most interesting result is the instability which the isotropic and homogeneous Universe 
acquires in the direction of the collapse toward the Big-Bang. For lo < 2/3 the density 
contrast diverges when approaching the cosmological singularity, i.e. for r\ — ► 0. In a 
non-viscous Universe the only perturbations which are able to generate the same kind of 
asymptotic instability are the tensor fluctuations (gravitational waves) whereas now the 
scalar perturbations destroy asymptotically the primordial Universe symmetry. The dy- 
namical implication of this issue is that an isotropic and homogeneous stage of the Universe 
can not be generated, from generic initial conditions, as far as the viscosity becomes smaller 
than the critical value associated to the condition lo = 2/3. 



VI. CONCLUDING REMARKS 



The main issue of our investigation is to have shown that the isotropic Universe acquires, 
backward in time, a regime of instability corresponding to sufficiently high values of the 
viscous parameter £o- Such a window of instability implies that, if the Universe was born 
sufficiently far from the homogeneous and isotropic stage, than the bulk viscosity (i.e. the 
absence of a stable thermal equilibrium) works against isotropization mechanisms and the 
inflation becomes the scenario from which a Robertson- Walker geometry arises (at least on 
a given scale). The explanation of this result is in the real physical meaning of the bulk 
viscosity: such viscous effects come out from the difficulty that microphysics finds to restore 
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the thermal equilibrium against the rapid Universe expansion. As a natural consequence of 
this physical context, bulk viscosity makes unfavored the establishment of an homogeneous 
stage from a general cosmological dynamics. On the other hand in a Robertson- Walker 
Universe, already settled down, we expect that, as we find, the viscous effects depress 
the density contrast because the particles inside the inhomogeneous fluctuations undergo 
dissipative processes which frozen the growth of the structures. Despite of the reliable 
feature of our result, the present investigation, as well as the whole previous literature on 
this subject, relies on a phenomenological ground; in fact the description of the viscous 
effect is based on the constitutive equation relating the viscosity coefficient to a power-law 
of the system energy density. This statement appears well-grounded, but nevertheless it 
requires to be carefully considered in a precise derivation of the viscosity coefficient from a 
real kinetic theory of matter [ijj ■ We will address for such a point in a further investigation, 
which will be aimed to yield an upgrading of the present cosmological issue. 
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